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ABSTRACT 

GRB1 10721 A was observed by the Fermi Gamma-ray Space Telescope using its two instruments the Large 
Area Telescope (LAT) and the Gamma-ray Burst Monitor (GBM). The burst consisted of one major emission 
episode which lasted for <~ 24.5 seconds (in the GBM) and had a peak flux of (5.7 ± 0.2) x 10~ 5 erg s _1 
cm~ 2 . The time-resolved emission spectrum is best modeled with a combination of a Band function and a 
blackbody spectrum. The peak energy of the Band component was initially 15 ± 2 MeV, which is the highest 
value ever detected in a GRB. This measurement was made possible by combining GBM/BGO data with LAT 
Low Energy Events to achieve continuous 10-100 MeV coverage. The peak energy later decreased as a power 
law in time with an index of —1.89 ± 0.10. The temperature of the blackbody component also decreased, 
starting from ~ 80 keV, and the decay showed a significant break after <~ 2 seconds. The spectrum provides 
strong constraints on the standard synchrotron model, indicating that alternative mechanisms may give rise to 
the emission at these energies. 

Subject headings: gamma-ray burst: general — gamma-ray burst: individual: GRB 11 0721 A — radiation 
mechanisms: thermal 
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1. INTRODUCTION 

Although the emission mechanisms active in the prompt 
phase of gamma-ray bursts (GRBs) are still under debate, 
there is much evidence that the photosphere of the relativistic 
outflow plays an important role in the formation of the ob- 
served spectrum (e.g. Lazzati & Begelman 2010, Ryde et al. 
2010, Guiriec et al. 2011, Vurm et el. 2011, Giannios 2011, 
Pe'er et al. 2012). Indeed, a strong contribution from the pho- 
tosphere was predicted on physical grounds in early works 
by Goodman (1986) and Paczyriski (1986) but this was not 
considered a viable model since the observed spectra are, in 
general, nonthermal. It was, however, realized that the photo- 
spheric emission should be accompanied by nonthermal emis- 
sion from optically thin regions (e.g. Meszaros et al. 2002), 
and that the photospheric emission can be enhanced and mod- 
ified from a Planck function by energy dissipation at moderate 
optical depths (Rees & Meszaros 2005). 

The Fermi Gamma-ray Space Telescope has made observa- 
tions that support this view. The very bright GRB090902B 
(Abdo et al. 2009) was observed to have a narrow and steep 
spectral component that may best be attributed to the pho- 
tosphere (Ryde et al. 2010, Zhang et al. 2011). At later 
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Fig. 1 .— Composite light curve of GRB110721A. (a-d) Count light 
curves from different energy ranges (Nal 6, 7, 9, 11 and BGO 1 
detectors), (e) Individual LAT photons. Time intervals are indicated 
by green lines; the red line shows the trigger time. Filled circles in 
(e) indicate > 90% probability of association with the GRB. The 
time is relative to the GBM trigger. 



times during the burst this spectral component widened into 
a broader feature. This implies that the spectrum of photo- 
spheric emission must be able to have a variety of shapes, 
not only a Planck function (Ryde et al. 2011). A probable 
explanation for this is the existence of sub-photospheric en- 
ergy dissipation. Furthermore, many bursts have shown signs 
of a subdominant photospheric component. For instance, 
GRB090820A exhibited two spectral peaks: a peak related 
to a blackbody spectrum with a temperature of ~ 40 keV, and 
a peak modeled by a Band function (Band et al. 1993) at 
E p ~ 1 MeV (Burgess et al. 2011). 

Only a few spectra having peak energies of a few MeV 
have previously been reported. Gonzalez et al. (2009) pre- 
sented GRBs observed over the energy range 0.02 - 200 MeV, 
through a joint analysis of data from BATSE and the EGRET 
TASC on the CGRO. The maximal value they found was 5 ± 2 
MeV in GRB981203. Similarly, the PHEBUS experiment on- 
board Granat, with an energy range of 0.1 - 100 MeV, iden- 
tified spectral breaks in other bursts at around 2.4 MeV in 
addition to the low energy break (Bar at et al. |2000"| ). Fermi 
observations of GRB090510 showed a peak energy of ~ 5 
MeV (Ackermann et al. 2010). 

In this paper, we study GRB110721A, which besides be- 
ing very bright during the first ~ 8 seconds (a few xlO -5 
erg s _1 cm~ 2 ) has an initial peak energy of a record break- 
ing 15 MeV - observations made possible by Fermi's excep- 
tional spectroscopic capability. Moreover, the 7-ray emission 
is dominated by a FRED (fast rise, exponential decay) pulse, 
which makes this burst ideal to compare with previous stud- 
ies of such pulses. Apart from the Band component, we iden- 
tify, with high significance, an additional component which 
we model using a Planck function. We interpret this as a pho- 
tospheric component. 

2. OBSERVATIONS 
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Fig. 2. — Spectral fit and residuals of the time integrated emission 
(-0.32-8.38 s) and the best fit model, Band + mBB. The top panel 
shows a vF u spectrum, the middle panel a count spectrum and the 
bottom panel the residuals of the fit. 

On 21 July 2011 the Gamma-Ray Burst Monitor (GBM; 
Meegan et al. 2009) and the Large Area Telescope (LAT; 
Atwood et al. 2009) onboard Fermi detected high energy 
emission from GRB110721A (GCNs 12187, 12188). The 
burst position was triangulated by the Interplanetary Network 
(IPN), which returned an error box centered in (RA, Dec) = 
332.46°, -38.63° (J2000) and approximately 1.3 x 0.4 deg 
wide (GCN 12195). The intersection of the IPN box with 
the LAT error circle gives an area of ~ 1200 square arcmin 
centered on (RA, Dec) = 333.2°, -38.5°. We adopted the 
latter position in our analysis. Figure [T] shows a composite 
count light curve from the various detectors on Fermi: The 
Nal detectors (8-900 keV), the BGO detector (200 keV-40 
MeV) and the LAT (P7V6_Transient class events, Atwood et 
al. 2009). The most energetic photon was detected at 4.50 s, 
had an energy of E — 6.3±0.6 GeV and was associated to the 
GRB with a high (> 0.9 using the gtsrcprob Fermi Sci- 
ence tool Q probability. The light curve of the burst is dom- 
inated by a single and exceptionally bright FRED-like emis- 
sion episode with Tgo = 24.5 s, the time during which 90 % 
of the emission is received. The peak of the energy flux occurs 
at 0.3 s relative to the GBM trigger with (5.7 ± 0.2) x 10~ 5 
erg s _1 cm -2 (over the energy range 8 keV - 1 GeV). 

In Fig. [T]we also present the light curve of the LAT Low En- 
ergy events (LLE, 30-130 MeV). The LLE data are produced 
from a non-standard LAT analysis which for bright sources 
provides large effective area at low energies, joining the LAT 
and GBM energy ranges (Pelassa 2011). An observation- 
specific response matrix is generated for each time range of 
data using a Monte Carlo simulation of the LAT, which in- 
creases the spectral capabilities of the Fermi data. We note 

1 http://fermi.gsfc.nasa.gov/ssc/data/analysis/ 



that the LLE light curve for this burst is peculiar since it peaks 
before, and its duration is significantly shorter than, the GBM 
light curve, in contrast to what is typically observed in other 
bursts {Fermi LAT Collab. in prep.). Finally, in Fig. [T]we ex- 
plicitly show the Nal light curve in the narrow energy range 
of 8 - 100 keV. Interestingly, at these energies the light curve 
differs in that a second peak appears at ~ 2 s. These facts will 
be discussed further in Sect. |4] 

We have performed a standard analysis with the RMfit 4.0 
package (Mallozzi et al. 2005) using TTE data from the Nal 
6, 7, 9, 11 and the BGO 1 detectors. We also include the LLE 
and LAT data. 

3. SPECTRAL BEHAVIOR 

A Band function poorly fits the integrated emission during 
the pulse (-0.32 - 8.38 s). The Castor C-statistic (C-stat; Dor- 
man et al. 2003) for the fit is 1078 for 618 degrees of freedom 
(dof). Based on the earlier detection of a blackbody compo- 
nent in other Fermi bursts (Ryde et al. 2010, Guiriec et al. 
2011, Burgess et al. 2011), we test whether the fit can be 
improved by including a blackbody. The fit improves signif- 
icantly to C-stat/dof = 901/616 (corresponding to > 5ct sig- 
nificance). Compared to the Band-only fit, the peak energy 
of the Band function shifts from E p = 1120 ± 60 keV to 
2400 ±170 keV. Since the blackbody component is expected 
to evolve, (see analysis below and Ryde & Pe'er 2009) its 
time-integrated spectrum is expected to be better character- 
ized by a multicolor blackbody (mBB; Pe'er & Ryde 2011). 
Using such a component instead of the blackbody yields a 
significant improvement of the fit C-stat/dof= 871/615. Our 
preferred fit for the time-integrated spectrum is thus the Band 
+ mBB model as shown in Fig. [2] 

We also analyze the time-resolved spectra using data up to 
130 MeV (i.e. the LLE data). In order to study the spectral 
evolution we need a temporal resolution as high as possible. 
However, since the signal-to-noise ratio (SNR) then also de- 
creases, there is a trade off in our choice of binning. We re- 
quire a SNR of 40 in the most strongly illuminated GBM de- 
tector (Nal 9) in the energy range 8.0- 100 keV. This ensures 
that the spectral fits are well constrained over the energy range 
where the blackbody component was found in the time inte- 
grated spectrum. This procedure divides the bursts into 8 time 
bins which are indicated in Table 1 . 

As for the time-integrated spectrum, we find that the ad- 
dition of a blackbody to the Band component improves the 
fit. In the first 7 bins the C-stat value decreases by a signifi- 
cant amount (see Table 1). Using an assumed Band function 
spectrum as the null hypothesis, and Band+blackbody as the 
alternative hypothesis, we determine the signicance of the ad- 
ditional blackbody component. The distribution of the test 
statistic AS was investigated through Monte Carlo simula- 
tions. We find that an improvement of AC-stat = 30 when 
adding a blackbody component to the fit, i.e., adding two de- 
grees of freedom, corresponds to a 10 -7 probability of not 
having a real blackbody component in the spectrum. Sim- 
ilarly, AC-stat = 20 corresponds to a 10~ 5 probability and 
AC-stat = 10 to a 10~ 3 (3ct) probability. In particular, the 
normalization of the blackbody component is constrained for 
the first 7 bins. We tested the robustness of this result by using 
another time bin selection based on Bayesian blocks (Scargle 
1998). The results are fully consistent. We therefore conclude 
that an extra component is significantly detected in the time- 
resolved spectra up to ~ 6 s. 

In Fig.[3]we plot the observed temperature of the blackbody 
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as a function of time. First, it is apparent that there is a strong 
evolution of the temperature. Second, there is a notable break 
in the decay of the temperature. The energy flux in the black- 
body component is approximately 5 % of the total flux and 
it peaks at around 2 s, when it reaches 10%. In the energy 
band around the peak it is much stronger, and well above the 
systematic errors ( 10%, Bissaldi et al. 2009). We note that 
this peak is notably different from the peak of the energy flux 
pulse, which occurs at 0.3 s, but coincides with the tempera- 
ture break. Moreover, it coincides with the second peak in the 
Nal count light curve (Fig. [T| which is mainly due to photons 
below 100 keV. This suggests that this second peak is associ- 
ated with the blackbody component. The temperature is well 
fitted with a broken power law (see Ryde 2004) in time with 
a break at 2.3 ± 0.2 s. Note that such a break does not appear 
in the decay of E p (Fig. 4). The power law indices before 
(after) the break are -0.30 ± 0.13 (-1.66 ± 0.15). This be- 
havior is quantitatively the same as what was found for the 
photospheric component in CGRO BATSE bursts (e.g. Ryde 
& Pe'er 2009). 

In order to study the evolution of the Band component we 
allow for a higher temporal resolution since it is the dominat- 
ing spectral component. We adjust the time bins to provide 
a minimum SNR of 30 in the counts spectrum from the Nal 
9 detector using the full energy range (8.0-860 keV), which 
increases the temporal resolution while still maintaining good 
spectral constraints. Since the burst is initially very strong in 
the LLE data and spectral evolution is initially very rapid, we 
divide the first time bin into two. This gives us 13 time bins. 

A striking feature in this burst is the unusually high value 
of the peak energy. The highest value is measured in the first 
time bin (-0.32-0.0 s). Here the peak energy is 15 ± 2 MeV, 
the low energy slope has a photon index a = —0.81 ± 0.08, 
the high energy power law has f3 — — 3.5+g'g. The fit has C- 
stat/dof of 679/608. Note that a blackbody component cannot 
be constrained by the data in this time bin, but we include it in 
the following bins. The peak energy of the Band component 
decreases monotonically with time; E p = A p \(t — t ) d , where 
d = -1.22 ± 0.13, and t = -0.46 ± 0.10 s (Fig. 0). At- 
tempting to fit the data before and after 2.3 s (the time where 
the temperature evolution shows a break) separately gives in- 
dices that are compatible within errors. We note that to is 
consistent with the onset of emission in the LLE range . 

After a few seconds E p reaches a few hundred keV, typical 
for other GRBs (Kaneko et al. 2006). Consistent results are 
obtained when using other temporal binnings. Figures [3] and 
[4] show the values found for both the coarser (red points) and 
finer (grey squares) time resolution. The two sets of points are 
clearly consistent. 

The binning used for the time-resolved analysis results in 
very few data above 100 MeV in each bin. We therefore also 
studied the time integrated spectra in three broad time inter- 
vals in order to increase the signal at the highest energies: 
[-0.325,1.5] s, [1.5,5] s, [5,20] s. For each of these bins 
an extrapolation of the best fit (Band+BB) model of the data 
below 100 MeV is consistent with the LAT data. The high- 
energy emission continues for more than 200 s after the burst 
trigger (Tgo for the LAT emission is 253 s) and decays as a 
power-law with index —0.95 ± 0.04, a typical value for the 
temporally extended emission observed in other bursts (Fermi 
LAT Collab. in prep.). 

4. DISCUSSION 
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Fig . 3 . — The blackbody temperature, kT, decays as a broken power 
law (fit function from Ryde 2004). Circles correspond to the binning 
based on the data below 100 keV, filled/open indicating 5a/3a sig- 
nificance of the blackbody component. Gray squares show results 
from the higher time resolution binning. 
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Fig. 4. — The evolution of E p of the Band function component. 
Symbols as in Fig. [3] 

The value of E p = 15 ± 1.7 MeV found in the initial time 
bin in GRB1 10721 A is the highest measured in a GRB spec- 
trum. Capturing such a high initial peak energy is a testament 
to the importance of acquiring good quality BGO and LLE 
data in bursts; this paper provides the first such realization of 
this contribution to the GRB paradigm. 

Bursts with hard spectra below their vF v peaks are fre- 
quently observed. Those with indices a > —1.5 below this 
peak cannot possess electrons that radiate synchrotron emis- 
sion in the expected fast cooling regime, within this spec- 
tral window; this is the so-called fast-cooling a index limit 
(Preece et al. 1998). GRB110721A is just such an exam- 
ple, and its spectrum can be consistent with optically-thin syn- 
chrotron emission only if the non-thermal electrons are cooled 
on timescales much longer than those on which they are in- 
jected/replenished. Demanding that the cooling break should 
lie well above 15 MeV, combined with the requirement that 
the observed vF u peak at 15 MeV in the first time window 
corresponds to the synchrotron peak, constrains the emission 
region to be smaller than typical photospheric radius scales 
(Ryde et al., in prep.), i.e. 10 10 — 10 13 cm when accounting for 
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TABLE 1 

Fit results using a Band function only and a Band+Blackbody model. 



Band function only Band+blackbody 
Time (s) -E pcak (keV) a g C-stat/dof fcT (keV) g peak (keV) a C-stat/dof 
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GRB bulk motions (e.g., see Rees & Meszaros 2005). This is 
a significant restriction on synchrotron emission models that 
invoke dissipation zones outside the photosphere. Therefore, 
it is desirable to entertain alternative explanations for the very 
"blue" emission seen in GRB 110721 A, e.g. Compton scatter- 
ing of thermal photons. 

The cleanest signature of the emission in GRBs is found 
in smooth FRED emission pulses (For d et al. |[T9 95). Ryde 
(2004) identified a few such pulses in the Compton Gamma 
Ray Observatory ( CGRO) BATSE catalogue which were con- 
sistent with having a Planck function spectrum throughout 
their durations over the observed spectral range 25-2000 
keV, suggesting a photospheric origin of the emission. It was 
further argued in Ryde (2005) that the spectral break in ap- 
parently nonthermal spectra can be interpreted as the photo- 
spheric emission peak. While the peak is modeled by a Planck 
function, an additional power law component is needed to ac- 
count for the nonthermal character of the spectrum. More 
importantly, the temperature of the blackbody in all pulses 
is observed to have a characteristic, recurring, behavior: it 
decreases following a broken power law in time (see further 
Ryde & Pe'er 2009). These results suggest that photospheric 
evolution in GRBs may exhibit well-defined characteristics. 

The evolution of kT seen here is similar to what was ob- 
served by BATSE over individual GRB pulses where a ther- 
mal emission component was identified (Ryde 2004). More- 
over, if this burst had been observed by BATSE it would 
be best modeled by a blackbody and a power law, which 
is apparent by limiting the energy range to that of BATSE. 
GRB 11 0721 A therefore confirms the BATSE results on the 
behavior of the photospheric emission. 

Finally, the two peculiarities of the 30-100 MeV emis- 
sion, its early onset and short duration, can be explained by 
the fact that the emission originates in the same component 
as detected by the GBM. This component initially dominates 
the LLE range due to the exceptionally high value of E p but 



moves into the GBM range as the burst evolves. Therefore the 
emission detected in the LLE range precedes, and has shorter 
duration than, the emission seen at lower energies; this mir- 
rors the behavior seen in BATSE bursts with FRED-like light 
curves (Norris et al. 1996). 

In summary, for GRB 11 0721 A we have shown that (i) ini- 
tially the peak energy has an extreme value E p = 15 MeV, 
the identification of which was made possible by the high 
quality LLE data. This value combined with the fact that 
a > —1.5 cannot be explained by the standard optically-thin 
synchrotron model, (ii) In addition to the Band function a nar- 
row spectral component (consistent with a blackbody) is sta- 
tistically required by the data, (iii) The spectral evolution of 
this blackbody is similar to what was found in BATSE pulses 
having a photospheric component, (iv) The blackbody com- 
ponent evolves differently compared to the Band component, 
(v) The second peak in the Nal count light curve coincides 
with the time when the blackbody component flux is the high- 
est. These facts provide strong evidence for the existence of 
photospheric emission in GRB 11 0721 A. 
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